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In  connection  with  the  first  goal;  l.e,  extending  the  basic  theory  to 


finite-element  problems,  one  report  has  been  written;  Block  Iterative 
Methods  for  Elliptic  Finite- Element  Equations.  This  69  page  report  is 
authored  by  Seymour  V.  Parter  and  Michael  Steuerwalt  (of  the  Los  Alamos 
National  Laboratory).  This  work  has  been  submitted  for  publication  in 
the  SIAM  Journal  on  Numerical  Analysis.  We  have  every  reason  to  believe 
it  will  be  accepted  for  Journal  publication.  This  report  contains  a 
general  theory  -  including  a  new  convergence  proof  for  a  large  class  of 
iterative  methods.  This  proof  does  not  require  that  the  matrices  involved 
be  symmetric  nor  does  it  require  that  the  matrices  be  of  poSitive-type  or 
inverse  monotone.  Of  course,  the  fact  that  these  matrices  arise  from  the 
discretization  of  elliptic  problems  is  essential.  The  second  part  of 
this  report  deals  with  the  application  of  this  theory  to  a  special  class 
of  problems. 

Let  {(x,y);  0<x,y<l}  be  the  unit  square.  Let 


L:=  l  a 


1j  dX^Xj 
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where 


c{x,y)  <  0  , 


be  a  uniformly  elliptic  (2nd  order)  operator  defined  on  fi.  Consider  the 
Dirichlet  problem 
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Lu  =  f  ,  ft 

u  =  0  ,  3ft  . 

Consider  also  a  finite-element  approach  for  the  numerical  solution  of  this 
problem.  Let  the  finite-element  method  be  based  on  the  finite-dimensional 
subspaces  given  by  the  tensor  products  of  Hermite  Cubic  Splines.  Let  us 
now  solve  these  equations  by  k-line  (horizontal  lines)  iterative  methods. 
The  results  are:  Let  Aq  be  the  smallest  eigenvalue  of  the  problem 

L$  +  Ac<J>  =  0  :  ft 
<{>  =  0  :  3ft  . 


Then 

A)  for  the  k-line  Jacobi  iterative  method 


p  *  1 
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B)  for  the  k-line  SOR  iterative  method  with  the  optimal  u> 


pb 


h  . 


The  method  is  quite  powerful  as  a  "comparer”  even  when  complete  results 
are  not  available.  For  example,  in  the  case  of  the  "point"  Gauss-Seidel 
method  we  can  not  give  a  complete  answer.  However,  we  have  the  following 
result.  Let 

q(x,y):=  ylf-  Ca11 (x,y)  +  a22(x,y)] 
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and  let  Aq  now  be  the  smallest  eigenvalue  of  the  problem 

Ld>  +  AQ<j>  =  0  :  n 

<J>  =  0  :  3ft  . 


Let  p  denote  the  spectral  radius  of  the  point  Gauss  Seidel  method. 
Then 

p  >  1  -  AQh2  +  0(h2)  . 


With  this  report  we  feel  that  this  topic  has  been  essentially  resolved. 
There  are  several  computational  projects  being  organized  in  connection  with 
this  work.  These  projects  will  highlight  some  of  the  salient  points  as  well 
as  Investigate  certain  delicate  questions  left  unresolved.  Professor  Parter 
will  report  on  this  work  as  part  of  his  "Retiring  Presidential  Address"  at 
the  SIAM  1983  Fall  meeting,  November  7-9,  1983  in  Norfolk,  VA.  Also,  he 
will  present  a  report  on  the  computational  results  and  related  topics  at 
the  Sixth  International  Symposium  on  Computing  Methods  in  Engineering  and 
Applied  Science  to  be  held  December  12-16,  1983  in  Paris,  France. 

In  connection  with  the  study  of  multi-grid  methods  a  first  basic  project 
has  begun.  There  Is  preliminary  theoretical  work  on  two  topics: 

(a)  1-dimensional  problems. 

(b)  MGR[0]  methods  for  problems  with  variable  coefficients. 

In  both  of  these  topics  the  theoretical  work  is  based  on  a  thorough  analysis 
of  "two-grid"  methods.  Computational  experiments  are  now  being  planned  and 
organized. 
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Finally,  on  the  topic  of  Conjugate  Gradient  methods,  a  small  experimental 
program  has  been  outlined.  This  program  is  based  on  certain  questions  which 
arise  in  the  analysis  of  "Block  Iterative  Methods  for  Elliptic  Finite-Element 
Equations".  In  simple  terms,  the  analysis  of  the  full  problem  was  reduced  to 
the  analysis  of  a  related,  but  simpler,  finite-difference  system.  The  conju¬ 
gate  gradient  program  will  study  these  two  problems  and  try  to  relate  the 
"rates  of  convergence"  for  the  conjugate  gradient  methods  in  these  cases. 


A 
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SUMMARY 


Reports  Completed: 

Block  Iterative  Methods  for  Elliptic  Finite-Element  Equations 
(with  M.  Steuerwalt) 


Work  in  Progress  (begun  during  the  period  6/15/82-6/14/83) 

Remarks  on  Multi -grid  for  1 -dimensional  Problems 
Remarks  on  MGR  Multi  grid  Methods 

Computational  Studies  on  Point  S.O.R.  for  Elliptic  Finite-Element 
Equations 

Computational  Studies  on  k-line  Iterative  Methods  for  Elliptic 
Finite-element  Equations 

Computational  Studies  on  k-line  Conjugate  Gradient  Methods  for 
Elliptic  Finite-element  Equations 


